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Report 

A.)  Statement  of  the  Problem: 


The  research  developed  mathematical  models  of  physical  phenomena  f^r  the  study 
of  complex  time-dependent,  nonlinear  partial  differential  equations  of  interest  in  nonlinear 
optics.  Applications  of  the  work  include  optical  communications,  optical  switching  and 
future  optical  computing. 

The  research  involved  the  development  of  mathematical  models  describing 
femtosecond  pulse  propagation  in  nonlinear  optical  media.  This  research  is  significant 
because,  as  the  duration  of  optical  pulses  shortens,  new  interactions  arise.  Recent 
advances  in  femtosecond  light  sources  makes  possible  the  study  of  new  phenomena. 
These  interactions  must  be  described  by  novel  nonlinear  partial  differential  equations.  A 
second  imponant  area  of  research  involved  the  interaction  of  multiple  beams  of  light  for 
intensity  dependent  ultrafast  optical  switching.  In  particular,  we  have  investigated  novel 
solitons. 


B.)  Summary  of  Results: 

Research  in  both  quantum  and  classical  areas  has  been  investigated.  We  have 
obtained  exact  solutions  to  coupled  higher-order  nonlinear  Schrodinger  equations.  This 
represents  the  first  work  in  this  significant  area.  These  equations  are  then  used  to  model 
femtosecond  all-optical  switching,  which  has  important  applications  in  the  optical 
computing  area. 
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Funhermorc,  using  the  time-dependent  Harytree  approximation,  we  have  derived 
the  first  investigation  of  quantum  effects  for  femtosecond  pulses.  This  new  work  may  lead 
to  a  greater  understanding  of  quantum  noise.  Our  results  describe  the  propagation  of 
femtosecond  solitons  in  nonlinear  optics.  These  solitons  travel  at  velocities  that  differ  from 
those  of  the  picosecond  solitons  obtained  from  the  standard  quantized  nonlinear 
Schrddinger  equation.  From  the  quantum  solutions,  we  find  that  the  soliton  experiences 
phase  spreading  and  self-squeezing  as  it  propagates. 
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Femtosecond  Solitons  in  Nonlinear  Optical  Fibers: 
Classical  and  Quantum  Effects 


F.  Singer,  M.  J.  Potasek,  J.  M.  Fang  and  M.  C.  Teich 

Departments  of  Physics,  Applied  Physics,  and  Electrical  Engineering 

Columbia  University 
New  York,  New  York  10027 


Abstract 

We  use  the  time -dependent  Hartree  approximadon  to  obtain  soludons  to  a  quandzed  higher-order 
nonlinear  Scirroedinger  equadon.  This  equation  describes  pulses  propagating  in  nonlinear 
optical  fibers  and,  under  certain  condidons,  has  femtosecond  soliton  soludons.  These  solitons 
travel  at  velocities  that  differ  from  those  of  the  picosecond  solitons  obtained  from  the  standard 
quantized  nonlinear  Schroedinger  equation.  Funhermore,  we  find  that  quadruple-clad  fibers  are 
required  for  the  propagadon  of  these  solitons,  unlike  the  solitons  of  the  standard  nonlinear 
Schroedinger  equation  which  can  propagate  in  graded-index  optical  fibers.  From  the  quantum 
solution,  SVC  find  that  the  soliton  experiences  phase-spreading  and  self-squeezing  as  it 


propagates. 


Soliton  solutions  to  coupled  higher-order  nonlinear  SchrOdinger 
equations 

R.  S.  Tasgal  and  M.  J.  Potasek 

Department  of  Applied  Physics,  Columbia  University,  New  York,  New  York  10027 

(Received  19  August  1991;  accepted  for  publication  30  October  1991)' 

A  set  of  coupled  higher-order  nonlinear  Schrodinger  equations,  which  describe  electromagnetic 
pulse  propagation  in  coupled  optical  waveguides,  is  formulated  in  terms  of  an  eigenvalue 
problem.  Using  that  result,  the  inverse  scattering  problem  is  solved  and  explicit  soliton 
solutions  are  found.  Additionally,  linear  coupling  terms  are  studied  systematically. 


I.  INTRODUCTION 

Optical  signal  processing  is  attracting  interest  be¬ 
cause  of  its  ultrafast  response  time.  Currently,  electro¬ 
optic  devices  generally  require  a  cumbersome  interface 
between  electronics  and  optics.  On  the  other  hand,  all 
optical  signal  processing  utili2ing  only  the  nonlinear  in¬ 
dex  of  refraction  results  in  partial  loss  of  the  signal  due  to 
nonuniform  intensity  distribution  within  the  pulse.  Soli- 
tons,  and  in  certain  cases  solitary  waves,  which  balance 
nonlinearity  and  dispersion,  can  improve  system  perfor¬ 
mance  due  to  their  remarkable  stability  properties. 

Slowly  varying  electromagnetic  waves  in  a  nonlinear 
medium  (an  optical  fiber,  for  example)  are  described  by 
the  nonlinear  Schrodinger  equation  (NLS).  In  order  to 
increase  bit  rates  it  is  necessary  to  decrease  the  pulse 
width.  As  pulse  lengths  become  comparable  to  the  wave¬ 
length,  however,  the  NLS  equation  becomes  inadequate, 
as  additional  terms  must  now  be  considered.  We  will  refer 
to  equations  which  include  these  terms  as  higher-order 
nonlinear  Schrodinger  equations  (HNLS). 

One  of  the  first  HNLS  equations  to  be  solved  exactly 
(by  Hirota'  in  1973,  two  years  after  the  simple  NLS 
equation  was  solved^)  and,  in  a  sense,  the  simplest  is 

iq,  -I-  i5q,  +  {P/2)q„  +  P\q\^q  -  ie(q,„  +  2fi\q\^q,)  =0, 

(1.1) 

where  fx  =  3  and  e  approaches  zero  when  the  pulse  width 
is  long  compared  to  the  wavelength. 

There  are  several  ways  to  generalize  the  HNLS  Eq. 
(l.l)  to  a  set  of  coupled  equations,  depending  on  the 
physical  situation  that  is  being  modeled.  A  fairly  general 
form  of  coupled  HNLS  equations  is 

iqi/  -f  /(5  +  -t-  5  _  )q\,  -t-  (P/l)qx„  -1-  ^(  |^i  |  ^ 

+  +  (^  +  -t-  A  _  )^i  -b  (A’^. 

+  iK _)qi-  ie[q\,„  +  n(\q\\^ 

+  y  kj  I  ^)^ir  +  +  rq*qi,)q\]=^^  ( i-2) 


klr  +  '(5+  —  5  +  (^/2)^2n  +  ^(ykl  1  ^ 

+  1^2  1^)92  +  (^+  ~  ^  -  )92  +  (^+  —  'A_  )^| 

-  [92m  +  M ( y  1 9i  I  ^  +  1 92 1  ^ )92:  +  /^ ( y9?9ir 

+  9*921)92!  =0-  (L3) 

A  nonlinear  directional  coupler  has  5_  =A^ 
=  A_  =  A_  =  0  and  K  +  9^.^  A  birefringent  single 
mode  fiber*  ’  and  rocking  fiber  rotator,’  in  which  the  fiber 
is  periodically  twisted,  have  A^  =  A_  =  0and5_,A_, 
A  ^  9^0,  where  y  is  a  function  of  the  ellipticity  angle  0 
and  two  material  parameters  a  and  b-. 


2o  +  2b  sin^  d 
la  -1-  6  cos^  0 


(Ref.7) 


in  optical  fibers  a=b. 

Equations  (1.2)  and  (1.3)  with  6_=A_ 

=  A+  =  A_  =  £  =  0  and  A  ^  9M3  describe  a  nonrelativ-  : 
istic  |x>son  field.'  In  a  weakly  relativistic  plasma,  nonlin¬ 
ear  coupling  of  two  polarized  transverse  waves  with  dis-  ■ 
persion  is  described  by  5+=5_=A+=A_ 

=  A^.  =  A_  =£  =  0.’  Also,  for  the  case 
5  _  =  )3  =  A+  =  A_  =  A+  =  A_  =  0  with  q\,qi  eR,  .  - 
Eqs.  (1.2)  and  (1.3)  are  a  pair  of  coupled  modified  ^ 
Korteweg-de  Vries  equations.  The  intennode  switching  :  'p 
term  A  _ ,  which  emerges  in  a  natural  way  from  the  math- 
ematical  derivation  below,  has  not  been  considered  in 
previous  soliton  work. 

Solitons  have  been  found  in  a  variety  of  (uncoupled) 
higher-order  NLS  equations.  Analytic  solutions  to  the 
simplest  NLS  equation — Eq.  (1.1)  in  the  limit  £— 0— 
were  discovered  in  1971  by  Zakharov  and  ShabaL^'® 
Hirota'  obtained  exact  soliton  solutions  to  the  HNLS  Eq. 

( 1.1 )  by  transforming  the  NLS  equation  into  a  homoge¬ 
neous  form  of  the  second  degree.  (While  this  approach 
produces  several  valuable  insights,  it  has  the  disadvan¬ 
tages  of  being  ad  hoc,  somewhat  hard  to  work  with,  and 
it  treats  the  higher-order  terms  and  NLS  terms  differ- 
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ently.)  Sasa  and  Satsuma  recently  discovered  soliton  so¬ 
lutions  to  a  more  complex  HNLS  equation."  The 
derivative"'"  and  mixed  derivative"'"  NLS  equations 
have  been  solved.  Painlevc  techniques  produce  other 
solutions."  Some  exact  soliton  solutions  may  be  found  by 
transformation  to  known  NLS  equations." 

Exact  solitary  waves  (which  may  or  may  not  be  soli- 
tons)  can  be  found  by  direct  substitution.'*  Approximate 
solitary  waves  can  be  found  by  various  techniques. 
Among  the  most  useful  are  numerical  computation"  and 
variational  methods.^  (Sometimes,  when  the  exact  solu¬ 
tions  are  not  known,  these  are  also  fruitful  approaches  to 
integrable  systems.^ 

Recently  the  coupled  NLS  equations,  without  higher 
order  terms,  have  been  the  focus  of  intense  attention." 
Manakov^*  found  explicit  soliton  solutions  to  the  coupled 
NLS  Eqs.  (1.2)  and  (1.3)  with  equal  self-  and  cross¬ 
phase  modulation,  but  without  either  higher-order  terms 
or  linear  coupling  5_=A^.=A_=  K+=K_ 

=  e  =  0,  y  =  1.  Elphick*  used  the  quantum  inverse  scat¬ 
tering  method  to  study  the  Manakov  equations  with  an 
added  symmetric  self  coupling  term  A  ^,60, 
5_— ^_=K_^=K_=£  =  0,  y=l.  Belanger  and 
Parc^’  found  a  transformation  that  reduces  a  set  of  cou¬ 
pled  NLS  equations  with  symmetric  linear  cross-coupling 
terms  K  +  5_=A^.=  A_=A'_=c  =  0,  y=l 

to  the  Manakov  case,  thereby  finding  solitons  with  peri¬ 
odic  energy  exchange  between  the  two  coupled  modes. 

Solitary  waves  in  the  coupled  NLS  equations,  includ¬ 
ing  noniiitegrable  (and  ccnsequently  soliton  destroying) 
terms  such  as  5_^,  and  absorption,  are  also  of 
interest.  We  mention  only  a  sampling  of  work  in  this 
field,  since  it  is  too  vast  for  a  thorough  survey  here. 
Christodoulides  and  Joseph"  discovered  exact  vector 
solitons  in  coupled  NLS  equations  with  a  birefringence 
term.  Pare  and  Flotjanczak"  found  analytic  solutions  us¬ 
ing  a  Lagrangian  variational  method.  Stability  analyses 
have  been  performed.^*'"  There  is  also  a  large  amount  of 
numerical  work.^" 

To  date  there  has  been  no  work  on  coupled  nonlinear 
Schrodinger  equations  with  higher-order  terms,  and  lin¬ 
ear  coupling  terms  have  not  been  studied  sys’''matically. 
In  this  paper,  using  the  method  of  Ablowitz,  Kaup,  New¬ 
ell,  and  Segur"  (AKNS),  we  formulate  the  coupled  NLS 
equations  in  a  more  systematic  way  than  has  been  done 
previously.  With  that  result,  the  coupled  NLS  equations 
are  generalized  in  a  very  natural  way  to  include  higher- 
order  terms  and  other  new  linear  coupling  terms.  The 
inverse  scattering  transform"'^*  is  then  straightforwardly 
carried  out,  yielding  explicit  solutions. 

II.  FORMULATION  OF  THE  EIGENVALUE  PROBLEM 

The  method  of  AKNS  is  begun  by  writing  the  not  yet 
fully  defined  eigenvalue  problem 


v,=  Tv, 

(-•P  9i  92  Y 

-9*  ip  oj,  (2.1) 

-9?  0  ip) 

v,=Zv, 

s 

Z,r  S  ^pp"-  (2.2) 

The  integrability  condition  for  Eqs.  (2.1)-(2.2)  is 

[T,Z]=0.  (2.3) 

Writing  each  of  the  nine  components  of  the  matrix 
explicitly  and  matching  terms  of  the  same  order  in  p 
yields  an  iterative  method  of  determining  Z: 

z{;' = (//2)(z{:V  "  -f  (z<;  ^ -  z<jj+ ")9, 


(2-4) 

Zi?>  =  (//2)(Z<;,^"  -f  (Zlr”  -Zjr '’)92 

-4r'’9i).  (2.5) 

Zi7'=  -  (i/inzil*  "  +  (Z<,7  +  "  -  Z<a"^  ")?• 

-Z<3^+"9?),  (2.6) 

-  (//2)(Z<7*"  -f  (Z<;+"  -z<3r”)9? 
-4rV).  (2.7) 

z\u=41^9^  +  (2.8) 

42"’=-47’9.-Or.  (2.9) 

(2.10) 

Z2j,’=-4;’92-Zi;’9f.  (2.11) 

^37*91 (2.12) 


and  also  four  equah'ties  in  the  zeroth  order  (1,2),  (1,3), 
(2,1),  and  (3,1)  matrix  components  of  the  integrability 
condition,  Eq.  (2.3): 

9u  -  -  (Zl?’  -  Z':S')9i  -b  45*92=0,  (2.13) 

92*  -  451  -  (45’  -  45’  )92  +  45*91 =o,  (lu) 
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-  Zl°V?  =  0.  (2.15) 

+  Z^°]  +  (^?’  -  -  Z\Vgr  =  0.  (2.16) 

Each  iteration  allows  five  constants  of  integration  in 
2^*’.  A  constant  times  the  identify  in  Z  or  7“  does  not 
affect  the  integrability  condition,  Eq.  (2.3).  There  remain 
four  possible  physical  degrees  of  freedom  for  each  term  in 
the  polynomial  Z. 

Setting  Z**’  =  0,  an  appropriate  choice  of  the  con* 
stants  of  integration  and  trace  yields 


/O  0  Ov 

Z<»=-8/c  0  1  0  , 

\o  0  .1/ 

/  0  qi  qi 

Z<J>=_4c  -qT  0  0 

^-q*  0  0 


-(f) 


0  0  0 
0  1  0 
0  0  1 


/l9il^+  l92l^ 

9i/ 

92/  \ 

1  /  Q\  i 

(  ° 

9i 

92\ 

1 

II 

-  l9il^ 

-?i«r|+2(f) 

1  -9? 

0 

\  9* 

“9192 

1 

0 

0/ 

-  ^i^f,  +  ^2/^2  -  ?iff  +  2(kil^+ 92«  +  2(|9i|^+ k2l^)92 


-?r«-2(kii'+  k2i')?r 

\-9?«-2(l?.|'+k2l')^? 


-  -  ?!?•) 
-  (?i/^2  -  9i9t) 


-  (ft<^r  -  ?29?,) 

-  (ft.^2  -  92^2^) 


l9il*+|92l*  9i,  92/  \  /O  0  °  \ 

9?,  -l9tl'  A^+A. 

92*,  -9(9?  -19217  \0  K^+iK. 


Wc  have  neglected  the  symmetric  group  velocity 
term  5  +  ,  which  enters  at  the  Z^‘ ’level,  since  it  can  be 
eliminated  by  a  trivial  change  of  variables.  Setting  Eq. 
(2.13)  equivalent  to  Eq.  (2.13)  and  Eq.  (2.14)  equiva* 
lent  to  (2.16)  forces  0,  A  +  ,  A_,  K_,  ceR. 

Insertion  of  Z^”’  into  Eqs.  (2.13)-(2.16)  yields  a  set  of 
coupled  HNLS  equations: 

'9ir+  (P/2)9i«  +  P(|9il^+  l92l^)9i  +  (A+  +  A,)^, 
4-  (Ar+  +/A:_)92-'f(9im  +  3(|9i|^+  |92l^)9i» 

+  3(^f9„  +  9?92/)9i]=0,  (2.21) 

'92r  +  (0/^)<}2ii  +  l9ll^  +  |92|^)92  +  (A  +  —  A_  )^2 
+  (K^  -iK_  )9i  -'f(92»/  +  3(|9iI*+ 

I-3(9?9i,  +  92  92/)92]=0.  (2.22) 

These  are  Eqs.  (1.2)  and  (1.3)  with  6^.=6_=0,  y 
=  1,  and  /i  =  3.  There  are  four  constants  of  integration 
introduced  in  Z*°’,  two  of  them  on  the  diagonal.  Equa¬ 


tions  (2.21)  and  (2.22)  therefore  contain  the  most  gen¬ 
eral  linear  coupling  terms  that  the  AKNS  formalism  al¬ 
lows  for  Eqs.  (2.l)-(2.3). 


III.  ELIMINATION  OF  THE  LINEAR  COUPLING 
TERMS 

Having  formulated  the  coupled  HNLS  equations 
above,  the  zeroth-order  constants  of  integration  may 
diagonalized  by  a  rotation: 

Ai>,=ArA-'Au,  (3.1) 

Ai»,=AZA-'Ai;,  (3.2) 

where 

,10  0  V 

A=  0  e-''^’^cos(fl/2)  e'^^’^sinffl/Q)  1 

Vo  -er<^’»sin(fl/2)  <“^’^c08(fl/2)/ 

(3.3) 
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tan(^)  =■ 


Un(e)  =  ^|K\  +  /:i/A  _ .  (3.5) 

This  is  equivalent  to  the  transformation 

9,=e<'^»»’[cos(d/2)9i  -sin(e/2)9il,  (3.6) 

9,=e-  <'^'«'[sin(6/2)9;  +  cos(6/2)9jl,  (3.7) 

with  primed  terms 

K\  =0.  (3.8) 


-.  =  ±  +  sign(A'_)=sign(A_). 

(3.9) 

and  all  the  other  terms  unchanged. 

The  linear  self  coupling  terms  A  ^  and  A'_  may  now 
be  removed  by  a  second  substitution: 


cannot  be  applied,  for  example,  to  a  periodically  twisted 
birefringent  fiber.  Neither  transformation  works  if  5  _  9^0 
or  Y^l. 


iV.  THE  INVERSE  SCATTERING  TRANSFORM 

Solitons  solutions  to  the  coupled  HNLS  equations 
may  now  be  found  using  the  inverse  scattering  transform. 

l9il>  1921-*^  ^  |r|-<aa  (which  implies  bright  soli- 
tons),  then  the  lost  functions  may  be  defined  as  the  eigen- 
functions  v  in  Eqs.  (2.1)  and  (2.2)  with  boundary  con¬ 
ditions 


—  [0|e  as  f—  —  00, 

\ol 


^  as  r-  -  00. 


(3.10) 


(3.11) 


^  as  r—  —  00, 

ili 


That  leaves  the  coupled  HNLS  equations  (omitting  the 
primes) 

'9u+  (P/2)9i„-f  P(|^,|^ -I-  |92p)9i-'fi9im 


^,1—  0  as  t—  00, 

lOi 


+  3(  1^1 1  ^  -f  19:1^)91,  +  3(9*9i,  -t-  9*92r)9i]— 0* 

(3.12) 

«9l:+  (^/2)92<r  +  P(|9ll^+  19:1^)92 -'c[92/fr 

+  3(19,1^  -I-  l92l^)92r  +  3(9?9„  +  9j92,)92]=0, 

(3.13) 

which,  in  the  limit  e->0,  is  the  Manakov  case. 

Belanger  and  Pare^  made  the  substitution 


9,=cos(A’  +  z)9',  -  /■  sin(Ar  +  z)95, 
92=  -  isin(A’^7)9j  -l-cos(if^7)92 


(3.14) 


to  eliminate  K  from  the  coupled  NLS  equations  (1.2) 
and  (1.3)  with  6_  =  A^  =  A_  =  K^  =c=0  and  y 
=  1.  In  contrast  to  the  transformation  given  by  Eqs. 
(3.6)-(3. 1 1 ),  Belanger  and  Fare’s  transformation  breaks 
down  (i.e.,  fails  to  eliminate  K  if  A  or 
although  it  does  work  for  the  other  linear  coupling  terms 
and  the  higher-order  terms.  Equations  (3.14)  and  (3.13) 


^,2—  1  as  r—  00, 

VOi 


as  /-«. 

For  peR,  T^=  —T.  Therefore 

^  ll^l(p.0^2(p.0]=0. 


(3.13)  where  is  a  solution  to  Eqs.  (2.1)-(2.3),  and 


=i\{p,t)^tj{p,t)  =5,y.  (4.8) 

The  lost  functions  (4.1)-(4.3)  are  related  to  the  Jost 
functions  (4.4)-(4.6)  by  the  scattering  matrix  a: 


i>uip,i)=  2  aij{p)i>,j{p,t), 
/-> 
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(4.10)  ^n^pJ)=ai  (p)^,i(p,0  +  au(p)\fr,2(p,0 


3 

X  a%(p)aji,(p)=8ij.  (4.11) 

k  =  \ 

Using  Eq.  (4.10),  aii(p)  ^  analytically  contin¬ 
ued  into  the  upper  half-plane  lm(p)>0;  and  a2j(p). 
a23(p),  ajzlp).  and  a^ip)  into  the  lower  half-plane 
Im(p)  <0.  From  this  and  the  unitarity  of  a  (4.11) 


a*  (p*)=det 


(ouip) 

U}2(p) 


a23(p)\ 

«33(p)  j 


(4.12) 


At  this  point  we  posit  that  On^p)  N  simple  ze¬ 
roes  at  the  points  p,,  pi,...,pff  in  the  upper  half-plane.  It 
will  be  shown  below  that  the  locations  of  the  zeroes  de¬ 
termine  (some  oO  the  physical  parameters  of  the  soli- 
tons. 

Introduce  an  integral  representation  of  the  lost  func¬ 
tions 


0,1  (P'U 


0,2(p.U 


(4.13) 


(4.14) 


(4.15) 


+  a23(p)0,3(p.^).  (4.19) 

0/3  (P.O  =a3i(p)0,i(p.O  +  a}2(p)0,2(p.f) 

+  a33(p)0,3(p.^)-  (4.20) 

Substitute  the  integral  representation  of  the  lost  func¬ 
tions  (4.13)-(4.I5)  into  Eqs.  (4.I8)-(4.20).  Operate  on 
Eq.  (4.18)  with 


where  goes  from  —  oo  to  oo,  over  p|...,p>  and  on 
Eqs.  (4.19)  and  (4.20)  with 


where  C_  goes  from  —  oo  to  oo,  under  p^,...,p%. 

This  gives  the  Gcl’fand,  Levitan,  and  Marchenko 
(GLM)  equations 

+  f*  +  s)f,2(t  +  s+r) 

Jo 

+  s  +  r)],  (4.21) 

0=Kt}\t,r)  +F2,(:  +  s)  +  jj  dslK^Ht.t  +  s) 


2) 


(3) 


The  functions  qy  and  qi  are  found  in  terms  of  K  by  sub¬ 
stitution  of  Eqs.  (4.14)  and  (4.15)  into  Eq.  (21): 


y-Fn{t  +  s  +  r)  -h  s)F2}(t  +  s  —  r)  ], 

(4.22) 

0=/rJj^(f,r) /'3i(f -hi) C  -h j) 

Jo 

XFjiit  -t-  s  +  r)  +  ArJ2*(t,t  -f  s)Fi2(,t  -f  j  —  r)  ], 


q,= 

(4.16) 

(4.23) 

q2= 

(4.17)  where 

To  find  and  we  first  recall  the  definition  of 
the  scattering  matrU  (4.9): 

0n(p.O=aii(p)0ri(p.O  +ai2(p)0r2(p.O 

■+-ai3(p)0,3(p.^).  (4.18) 


aiy(p) 

aii(p) 


(n. 
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dse^^ 


s 


=  -il  +  ^ 


a.Xp) 


<(P*) 


(-0. 


to  give  ?  set  of  iN  linear  equations  in  1  +  —  p* ), 

Kriit,p„),  and  k  r3(t,p„),  with  m  = 


1=2,3. 

Note:  the  z  dependence,  which  was  previously  implicit  in 
a,  is  now  written  explicitly  in  C,y. 

The  residues’  contribution  to  Fjj  corresponds  to  soli- 
tons;  the  Fourier  transform  part  corresponds  to  radiation. 
Since  we  wish  to  obtain  soliton  solutions,  we  will  neglect 
the  latter.  Performing  the  integrations,  Eqs.  (4.21)- 
(4.23)  become 

s 

0=Kll\t,r)-i  1  e^"<'^'>IC,2(z,p,)^,,(r,p,) 

1 

+  C,3(r,p,)^,3('^Jl.  (4.24) 


1  =  1  +^„(f,  -pj)  -  I  -i—T[Q2(z4>.)KnU4>.) 
ffsl  Pm  Pr 

+  C,3(r,pJ^,3(/,pJ],  (4.29) 

^  V 

0  =  ^r2(^P/n)  +  X  ~  ^^2l(z^*) 

n=\  Pm~  Ph 

X(1+^^,(/,-P:)),  (4.30) 

o=fCfi(t,p„)  +  X  “ 

«=1  pm  —  pH 


S 

Q=K^,^\t,r)-i  1  e-^*{t  +  r)C  (24>*) 

X(H-^„(A-P?)).  (4-25) 

0=A;3'>(r.r)-/  X  e-«’^+'>C3,(r,pJ) 

X(1 +^„(^-p?)).  (4-26) 

where  a  hat  denotes  the  operator 

;4(r,p)=  r“*e-^-4<'>(^r  +  J). 

Jo 

To  find  g,  and  qi,  set  r=t  in  Eqs.  (4.25)  and  (4.26); 

q,= 

N 

=  -2/  X  €-'^*‘C^^{z4>:)(\+k,,U,-p*)\ 

R^\ 

(4.27) 

N 

^~2il  €-'^?'C„(r,p*)(l  +^„(r.-p;)). 

(4.28) 


X(l+^„(^-p;))-  (4.31) 

The  z  dependence  of  the  Jost  functions  (4.1)-(4.6) 
may  be  determined  from  Eq.  (2.2).  For  simplicity,  we  use 
the  z  dependence  of  the  Jost  functions  in  the  limit  1 1|  — 0 
to  determine  the  z  dependence  of  and,  consequently, 
Cjj  (for  all  r): 


cc  exp{/(j3/2)  (I  +  iTj)*  -E(^  +  ii?)^]^} 


oc  exp{77[  (^/2) (2|)  -  c(3f*  -  7/^) ]r} 

Xexp{/[  (^/2)(^^  -  V)]r 

(4.32) 

Substitution  of  the  z  dependence  of  the  C,yS,  Eq. 
(4.32)  into  Eqs.  (4.29)-(4.31 ),  and  substitution,  in  turn, 
of  the  appropriate  results  of  Eqs.  (4.29)-(4.3 1 )  into  Eqs. 
(4.27)-(4.28)  gives  exact  soliton  solutions. 

V.  THE  ONE-SOLITON  SOLUTION 


Now  return  to  Eqs.  (4.24)-(4.26).  Substitute  r=f 
+  5  and  operate  with 


The  simplest  nontrivial  soliton,  found  by  setting  N 
=  1  above,  is  a  single  solitary  wave: 
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1  -  .fC,3(p)C3, (/.*)]  ’ 

1  -  [e^‘^>’-P’^‘/{p-p*)^][Cnip)C2dp*)  +  Cu(p)Cii{p*)]  ' 


(5.1) 

(5.2) 


On  substitution  on  the  z  dependence  Eq.  (4.32)  and 
2p  =  ^  +  ir],  and  some  algebra,  we  may  express  and  ^2 
in  the  form 


qi=sm{a)e">’-q", 

(5.3) 

q2=cos{a)e~‘*-q’. 

(5.4) 

where 

^'=7)sech{77{r-ro-f-  [(^/2)(2|’)  -£(3|’^-7/^)]z}) 

Xexp(  -/{|■(^-ro)  +  [(^/2)(|’^-t)^) 

*  / 

-W)\2 +  <!>+}} 

(5.5) 

tan(a)  =  |  Ci2(0,p)  |/1  CnCO,/?)  [ . 

(5.6) 

Finally,  the  transformation  given  by  Eqs. 
(3.9)  gives  the  unprimed  q^  and  q^: 

(3.6)- 

[cos  ( 0/2 )  sin  (a)  ^  ^  ’ 

—  sin(0/2)cos(a)e“^^'^-*'^^-  ^]q. 

(5.7) 

^2=e“^'^^’’’[sin(0/2)sin(a)e'^-*'''^-’ 

-f  cos(0/2)cos(a)e~'^'^-^'''^-^]9, 

(5.8) 

where 

^=7)  sech[7)0  -  fo  +  [ ((3/2)  (2^)  -  £(3^^  -  ]z}) 

Xexp[  —  /{^(t  —  to)  +  [  -  ^  +  ((3/2) (1^^  - 

•17^) 

-£^(^^-37)^)]z-f-i^  +  }j; 

(5.9) 

a,  (f>^,  (f>_,  and  Iq  are  free  (real)  parameters,  as  are  the 
components  of  the  eigenvalue  ^/2  and  q/2.  Recall  that 
tan(^)  =K_/K^, 

taniO)  =ylK\  +.^i/A_, 
and 

A'_  =  ±  +K\  +k'‘-_. 


The  other  parameters  are  given  in  the  coupled  HNLS 
equations. 

In  conclusion,  we  have  obtained  bright  soliton  solu¬ 
tions  to  a  generalized  set  of  coupled  higher  order  nonlin¬ 
ear  Schrodinger  equations.  Higher-order  and  NLS  terms 
are  treated  the  same  way.  Also,  we  found  a  transforma¬ 
tion  that  eliminates  all  of  the  four  linear  coupling  terms 
that  the  AKNS  formalism  allows  for  this  problem.  Fu¬ 
ture  papers  will  focus  on  dark  soliton  solutions  and  mixed 
dark  and  bright  soliton  solutions. 
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.Abstract.  We  use  the  time-dependent  Hartree  appro.ximation  to  obtain  the  solution  to  the 
quantum  higher-order  non-linear  Schrodinger  equation.  This  equation  describes  femto¬ 
second  pulses  propagating  in  non-linear  optical  fibres  and  can  have  soliton  solutions. 
These  solitons  travel  at  velocities  that  differ  from  the  picosecond  solitons  obtained  from 
the  standard  quantum  non-linear  Schrodinger  equation.  We  find  that  these  femtosecond 
solitons  cannot  propagate  in  graded-inde.x  fibres;  rather,  they  require  quadruple-clad 
fibres.  This  is  the  first  investigation  of  quantum  effects  in  femtosecond  solitons  to  our 
knowledge. 


There  is  considerable  interest  in  the  non-linear  Schrodinger  equation  (nls)  in  terms  of 
both  classical  and  quantum  phenomena  (1-6).  In  particular  it  has  been  used  e.xtensi- 
vely  to  model  the  propagation  of  pulses  in  non-linear  optical  fibres;  however,  the  nls 
is  generally  not  valid  for  pulses  with  durations  shorter  than  the  picosecond  time  scale/. 
Yet  the  recent  development  of  optical  sources  that  generate  pulses  in  the  femtosecond 
domain  makes  possible  the  e.xploration  of  many  new  phenomena.  Therefore  the 
investigation  of  solitons  arising  from  the  higher-order  nls  (hnls),  which  can  be  used 
in  the  femtosecond  time  domain,  is  of  interest. 

One  of  the  simpletst  hnls  is  (7j 


da  d'(p  d^<p  d(p 

i-  +  ^  +  2C|0|-9  +  W^  +  ipCl<pP-  =  O. 


(1) 


where  C,  d  and  p  are  constants.  We  follow  the  conventional  notation  in  the 
mathematical  literature,  which  uses  f  and  .r  to  represent  normalized  space  and  time, 
respectively.  This  equation  gives  rise  to  soliton  solutions  when  p  =  6c/  [7].  Equation  (1) 
reduces  to  the  nls  for  p  =  c/  =  0. 

In  certain  circumstances  the  hnls  can  be  used  to  describe  femtosecond  pulses 
propagating  in  optical  fibers;  these  are  outlined  in  [8]  and  described  in  detail  by  us  in 
[9].  Using  experimental  fibre  parameters  to  evaluate  the  physical  parameters  in 
equation  ( 1)  we  find  that  the  pulse  width  must  be  below  200  fs  for  wavelengths  in  the 
1.4S-1.57ccm  region  in  order  for  d  and  p  to  become  significant.  In  addition,  the 
dispersion  parameters.  and  Py,  given  by  the  second  and  third  derivatives  of  the 
propagation  constant  respectively,  evaluated  at  the  carrier  frequency  co,).  must  be 
nesative.  This  necessitates  a  quadruple-clad  fibre  rather  than  the  typical  graded-inde.x 
fibres  used  in  calculations  and  experiments  to  date.  This  is  a  significant  feature  of  our 
results  [9].  The  soliton  self-frequency  shift  (ssfs)  [10. 11]  may  be  an  important  effect 
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when  considering  femtosecond  solitons.  However,  we  use  the  numerical-beam  propa¬ 
gation  method  to  show  that  at  distances  required  for  the  quantum  effect  to  be 
observed  the  effect  of  the  ssfs  on  the  soliton  described  by  equation  (I)  can  be 
neglected  [9]. 

In  the  case  of  optical  soliton^^  represents  the  nornj|i(ized  envelope  of  the  electro¬ 
magnetic  licld.  The  quantities  C  and  d  arc  given  by 


- 


n^ojaa'I' 


3a|^:r 


and  p  is  a  parameter  involving  the  frequency-dependent  index  of  refraction  and  the 
frequency-dependent  radius  of  the  mode  of  the  fibre  [8].  n.  is  the  non-linear  index  of 
refraction,  a  is  the  w-idth  of  the  pulse  rln.-arign.  /  is  the  peak  amplitude  of  the 
pulse  and  c  is  the  speed  of  light. 

The  general  solution  of  equation  (1)  has  the  form  (7j 

(p  =  ip„  sech[f  (.r  -.ru)  -r^r|  exp{i[ y{x  -.to)  (3f]}.  (3) 

where  r.  y  and  b  are  constants  and  is  the  zero  of  time.  Substituting  this  in 
equation  (1)  yields  the  following  relations 

|<p„l’  =  £^'C  d  =  £'-y'"2dy£--rdy^  ^  =  £{2y-r  d£' -2dy').  (4) 

We  proceed  by  considering  the  quantum  version  of  equation  (1)  from  a  mathe¬ 
matical  point  of  view.  In  (9|  we  examine  the  physical  aspects  of  this  problem  in  detail 
and  describe  the  role  played  by  other  effects  such  as  the  ssfs.  The  initial  portion  of  our 
analysis  closely  follows  that  of  Lai  and  Haus  (5|  for  the  .vis.  To  obtain  the  quantum 
version  of  equation  (1).  the  quantities  0(r. .t)  and  tp'it.x)  are  replaced  by  the  field 
operators  <p(r,.r)  and  (p'(t.x),  which  satisfy  the  boso^  commutation  relations 

[(?(;.  .r').  9*(r.  .r)J  =  <)(.r-.r')  (^(r..r').  (p{r.x)|  =  x').  .r)]  =  0  (5) 

where  <p(.i.x)  and  (p*(r.  x)  are  the  photon  annihilation  and  creation  operators, 
respectively,  at  t  and  x. 

The  quantized  equation  can  be  written  as 

i/i^<p(r.x)  =  [(?(r..r)./?I  (6) 


1/e 

f 


^  =  J (pr(r..r).p,(f..r)dLr-C J (p‘{t.x)(p'it,x)<p{t.x)ip(i,x)ix 


(7) 
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We  solve  equation  (13)  using  the  time-dependent  Hartree  approximation  [12],  We 

inr*  1  Hnrtrpf*  \vn vpfiini'rinn 


define  a  Hartree  wavefunction 

. 

where  <I>,  has  the  normalization 
|<t»,(.t,r)|-dt  =  1. 


/-I 


/' 


(15) 


The  functions  O,  are  determined  by  minimizing  the  functional 


. 

E  ^(•f,-  r.)^2CT6iCd^^  /;”’{.r . .  r)  dr,  •  •  •  dr,.  (16) 

which  provides 

i  — V  2C(/i  -  1)  j  4)„i-4),  +  i^-^"  ^  6iCd(/i  -  1)14>,P  — '  =  0.  (17) 

at  ax  — o.r'’  —  o.r 

This  is  identical  to  the  classical  knls  given  in  equation  (1),  with  C  replaced  by 
C{n  -  1).  Thus  the  solution  to  the  quantized  femtosecond  soliton  equation  is  obtained 
directly  from  equations  (3)  and  (4): 

^..(.r.  t)  =  [C(n  -  !)]■"■ «  sech{«((.t  -  .r„)  -r  (-3dy*  +  rff*  +  2>')r)} 

X  exp[  -  \{dy'  -r  3(/y£-  +  y*  -  £  *)/  -f  iy(.r  -.r.j)].  ( IS) 

The  normalization  condition,  equation  (11).  gives 

f  =  i(/i-l)C.  (19) 

Substituting  equation  (19)  into  (IS)  leads  to 

d',„(.t.  i)  =  :('>  -  l)'''C''-  sech{:(«  -  l)C((.r  -.r,,)  +  ( - 3dy-  +  id(n  -  1)-C-  -r  2y)f]} 

X  exp{(id7'^-  ^y(/i  -  1)'C--  iy*-!-  -  I)’C-]r  +  iy(.r-.r,,)}.  (20) 

The  Hartree  product  eigenstates  are.  using  equations  (10)  and  (14), 


\n.y.t)H  = 


Vn'. 


<I>„(.r.f)<p*(.r)dr  |0> 


:  |l 


(21) 


A  superposition  of  these  states,  using  a  Poissonian  distribution  of  n  for  a  coherent- 
state  pulse,  gives 


(22) 


I 


where  |a,l*  =  «.,  is  the  mean  photon  number. 
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Figure  I.  Plots  of  ihe  qussiproftability  density  Qla.x.i)  against  the  real  and  Imaginary 
parts  of  Cl  for  <3a=  -1.  C  =  0.25.  i/  =  0.25.  y^O.  =  0  and  (u)  /«0  and  (6)  r  =  0.l. 


The  quasiprobability  density  for  the  amplitude  of  the  envelope  of  the  field  is  defined 
as 

Q(a..r.0  =  l(a..r|->O|’  (23) 

wher. 


'  a’ 


la.. r) a e ^ 


is  a  local  coherent  state  at  the  time.r.  Substituting  equation  (22),  with  (20)  and  (24), 
into  (23)  gives 

2(a..r.r)=c''""'‘"'V 


X  scch 


f/i-I 


(.r-.r„)+  -Mly^d - ; — C-  +  2>'  r 


:e.xp 


i^ndy  '‘  -  i/t(n  -  1)*  \dyC  '  - 


iny- in  ■ 


CMf 


I 

I 

I 

I 

I 


■r  i’/ri(.t  -  .r„) 


(23) 


In  figure  1  we  illustrate  how  this  quantity  changes  as  the  soliton  propagates  in  space. 
We  have  ignored  the  n  dependence  of  the  amplitude  and  kept  it  in  the  phase.  We 
observe  phase  spreading  similar  to  that  in  the  .nls  case  (5, 6j. 
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